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In this paper we develop the spectral theory of the fractional 
Brownian motion (fBm) using the ideas of Krein’s work on continu¬ 
ous analogous of orthogonal polynomials on the unit circle. We ex¬ 
hibit the functions which are orthogonal with respect to the spectral 
measure of the fBm and obtain an explicit reproducing kernel in the 
frequency domain. We use these results to derive an extension of the 
classical Paley-Wiener expansion of the ordinary Brownian motion 
to the fractional case. 

1. Introduction. Let X = (Xt)t>o be a fractional Brownian motion (fBm) 
with Hurst index H G (0,1), that is, a continuous, centered Gaussian process 
with covariance function 

EXsXt = + t‘^^ -\s- t\‘^^). 

Say the process is defined on the probability space P), and for some 

fixed time horizon T >0, define the linear space TSt as the closure in (P) of 
the (complex) linear span of a collection of random variables {Xt :t € [0, T]}. 

So-called linear problems for the fBm are problems in which it is re¬ 
quired to find elements of the Hilbert space TCt with certain specific prop¬ 
erties. In the 1960’s the basic linear problems like prediction, interpola¬ 
tion, moving average representation etc. were treated by various authors; 
see Molchan (2003) for an overview of these contributions. However, these 
results did not become widely known. Many of them were rediscovered 
during the last decade when new application areas like telecommunica¬ 
tion networks and mathematical finance stimulated a renewed interest in 
the fBm. Recent contributions dealing with linear problems can be found, 
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for instance, in Gripenberg and Norros (1996), Decreusefond and Ustiinel 
(1999), Norros, Valkeila and Virtamo (1999), Nuzman and Poor (2000) and 
Pipiras and Taqqu (2001). 

There exist several representations of the fBm that give insight into the 
structure of the linear space TIt- An important example is the spectral 
representation 

r (giAt _ iVg-iAs _ 

]EA,At= / ^- ’-f,{dX) = {eues)f,, 

JR A 

where fi{dX) = (27r)“^ sin(7r//)r(l + 2H)\X\^~‘^^ dX is the spectral measure 
of the fBm and et{X) = (e*^* — l)/iX [see, e.g., Yaglom (1987), page 407 or 
Samorodnitsky and Taqqu (1994), page 328]. If we define Ct as the clo¬ 
sure in T^(/r) of the (complex) linear span of the collection of functions 
{et'.t G [0,T]}, this representation gives rise to an isometry between TCt and 
the function space Ct, determined by the relation Xt < —e*. We can use 
this spectral isometry to reformulate a linear problem for the fBm in spectral 
terms. It then becomes a linear problem in the function space Ct, which has 
the advantage that we have mathematical tools like Fourier-type techniques, 
at our disposal. 

In this paper we present new results regarding the fine analytical structure 
of the frequency domain Ct- In particular, we exhibit certain “orthogonal 
functions” with respect to the spectral measure // of the fBm and we obtain 
an explicit reproducing kernel for Ct, turning it into a reproducing kernel 
Hilbert space (RKHS). To illustrate the significance of these new frequency 
domain results for the fBm, we apply them to derive a generalization to 
the fractional case of the classical Paley-Wiener expansion of the ordinary 
Brownian motion [cf. Paley and Wiener (1934)]. In spectral terms, obtaining 
a series expansion translates to finding an orthonormal basis of the space Ct- 
We achieve this by using the RKHS structure and the explicit expression 
that we have for the reproducing kernel. 

It is well known that orthogonal polynomials on the unit circle are very 
useful in the spectral analysis of stationary time series. They can be used to 
solve problems like prediction and interpolation, and are also useful in con¬ 
nection with likelihood estimation and testing [see, e.g., Grenander and Szego 
(1958)]. In a classical paper, Krein (1955) introduced certain continuous ana¬ 
logues of orthogonal polynomials on the unit circle. We refer to Akhiezer and Rybalko 
(1968) for a more elaborate treatment, including detailed proofs of Krein’s 
statements. As was shown by Kailath, Vieira and Morf (1978), Krein’s or¬ 
thogonal functions play the same role in the spectral theory of continuous¬ 
time processes with stationary increments as the orthogonal polynomials 
do in times series theory. For a certain class of processes with stationary 
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increments, Kailath, Vieira and Morf (1978) pointed out how the orthogo¬ 
nal functions describe the structure of the frequency domain. In the present 
paper we develop the spectral theory of the fBm along the same lines. 

The results of Akhiezer and Rybalko (1968) and Kailath, Vieira and Morf 
(1978) highly depend on the “signal plus white noise” structure of the pro¬ 
cesses that they consider. It turns out, however, that Krein’s ideas are also 
applicable for the fBm, which is not of the latter type. The key point is 
that the fBm can be “whitened” in the sense that integration of an ap¬ 
propriate deterministic kernel with respect to the fBm yields a contin¬ 
uous Gaussian martingale, the so-called “fundamental martingale.” This 
was hrst proved in the 1960’s by Molchan [cf. Molchan (2003)]. For alter¬ 
native, more recent approaches see also Decreusefond and Ustiinel (1999), 

Norros, Valkeila and Virtamo (1999), Nuzman and Poor (2000), Pipiras and Taqqu 
(2001) or Dzhaparidze and Ferreira (2002). Conversely, it is well known that 
by integrating a certain deterministic kernel with respect to the fundamental 
martingale, we can recover the fBm. Since these results play an important 
role in this paper, their precise statements will be recalled in Section 2. 

The whitening and moving average formulas for the fBm provide us with 
the starting point for the development of the spectral theory. They give rise 
to a Hilbert space isometry U between Ct and the space T^([0, T],V), where 
V is the variance function of the fundamental martingale. In the ordinary 
Brownian case H = 1/2 this isometry is simply the Fourier transform. We 
will show that, for H 7 ^ 1/2 it is also a Fourier-type integral transformation, 
and obtain an explicit expression for the Fourier kernel in terms of Bessel 
functions. Using this Fourier kernel, we then introduce a function St that 
will turn out to be a reproducing kernel on Ct- An explicit expression for 
this kernel will be derived in a number of steps. First we shall use the Bessel 
differential equation to prove that the properly normalized Fourier kernels 
satisfy Krein’s continuous version of the recurrence relation for orthogonal 
polynomials. It will then be rather straightforward to obtain a Christoffel- 
Darboux-type formula for St- In combination with the expression for the 
Fourier kernel, this will lead to an explicit formula for the reproducing kernel 
on Ct- This program is carried out in Sections 3-6. 

In Sections 7 and 8 we use the new results on the structure of the space Ct 
to derive an extension to the fractional case of the classical Paley-Wiener 
expansion of the ordinary Brownian motion. We will first use the reproducing 
kernel to find a suitable orthogonal basis of Ct- By transporting of this basis 
to the space TCt, we will prove that the fBm admits the series expansion 


E 


2iuJn 


1 

— Zn, 


fE[0,l], 


where the ujn are the real-valued zeros of the Bessel function Ji-h and the 
Zn are independent, complex-valued Gaussian random variables with zero 
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mean and a variance that can be expressed explicitly in terms of Bessel 
functions and their real zeros. Using the fact that = vV^sinz, 

it can be seen that for H = 1/2 this indeed reduces to 


E 


g2m7rt _ 

2imr 


-Zn 


tG[0,l], 


with the Zn i.i.d., standard Gaussian. This is the expression that Paley and Wiener 
(1934) used as the definition of the standard Brownian motion. 

We will also briefly consider questions like the rate of convergence of the 
Paley-Wiener expansion, and possible extensions to the fractional Brownian 
sheet. In particular, we will argue that the expansion is rate-optimal, in the 
sense of Kiihn and Linde (2002). This is obviously a desirable feature if the 
expansion is used for simulation purposes and is also relevant in connection 
with the small ball problem for the fractional Brownian sheet [see, e.g., 

Li and Linde (1999) and Li and Shao (2001) for the precise connections]. 


2. Auxiliary facts and notation. The spectral representation can be used 
to define a stochastic integral with respect to A of a large class of deter¬ 
ministic integrands. In this paper we denote the indicator function l(o,t) of 
the interval (0,t) simply by If Using this notation, as well as the previous 
notation et(A) = (exp(fAt) — l)/fA, we may write e* = If Here and elsewhere 
below we adopt the usual convention to denote the Fourier transform of a 
function / G L^(M) by /, that is, 

/(A)= / f{x)e^^^dx. 

Jr 


Now consider the class of functions It = {f & L^[0,T]:/ G L^(M,B(M),//)} 
and endow it with the inner product Then the spectral 

representation can be written as KXgXt = (Is) particular, the map¬ 

ping If —> Xf extends to a linear map I :It ^ 'Hr with the property that 

m!) = X,andt„r/.<,eXT. 


We denote the random variable /(/) by / f dX or Jq f{t) dXt, and call it 
the integral of / with respect to X. We note that, in general, not every 
element of TCt can be represented as such an integral since for H > 1/2 the 
space It is not complete [see Pipiras and Taqqu (2001)]. 

Let us now introduce an integral with respect to X which plays an im¬ 
portant role in this paper. For t>0, define the kernel mt by 


( 2 . 1 ) mt{u) 


_ t _ 

2HT{H + l/2)T{3/2- H) 
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where T denotes Euler’s gamma function. Then, for every t £ [0,r], it holds 
that mt €It, and in view of Poisson’s integral formula for the Bessel func¬ 
tion [e.g., Watson (1944), Section 3.3], it is not hard to see that the Fourier 
transform tht of mt is given by 


(2.2) mt{X) 


2Err(Er +1/2) 



_ ^ _ 

mvifl + i/2)r(2 - 


^2-2H^ 


A 7^0, 


A = 0, 


where Ji-h is the Bessel function of the first kind of order 1 — H [for details 
see Dzhaparidze and Ferreira (2002), Proposition 2.2]. On evaluating rht at 
A = 0, one has to take into consideration the basic property 


(2.3) 




1 

2^r{u + i) 


as z - 


0 


of the Bessel function. For convenience, we introduce a special notation 
for the constant that occurs in the second line of (2.2). In the literature this 
constant is often given in an alternative form, namely. 


(2.4) 


2 r{3/2-H) 

^ 2HT{H + l/2)T{3-2H)' 


The identity of the two expressions is a result of Legendre’s duplication 
formula for the gamma function. 

Now, for t £ [0,T], we can consider the random variable 


(2.5) ^ ^ fnt{u) dXu 

in TCt- As is proved in Dzhaparidze and Ferreira (2002), Theorem 2.3, it 
holds that KM^Mt = {ms,mt)ti = ^sAt(O). This shows that the process M 
defined by (2.5) is a continuous Gaussian martingale with bracket (M) = 
m(0). For convenience, this variance function will be denoted by V, so that 
due to (2.2), we have 

( 2 . 6 ) Vt = EM^ = 


for all t > 0. Following Norros, Valkeila and Virtamo (1999), we call the 
process M the fundamental martingale. 

Next we recall the moving average representation of the fBm, which is the 
converse of (2.5). Let xt be defined by 

xt{u) = - Ijt - v)^-^/^dv^-^/^yt{u). 

Then it holds that 

(2.7) W= [\t{u)dMu 

Jo 
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for all t>0, where M is the fundamental martingale. More precisely, the 
process on the right-hand side defines an fBm with Hurst index H, so, in 
particular, we have that xt G L^([0,r],H[0,r], H) for all t G [0,T] and 

rsAt 

(2.8) EXsXt= Xs{u)xt{u)dVu = {xs,xt)v 

Jo 

for G [0,r]. We therefore consider the space /Ct, defined as the closure 
in i/^([0, T],B[0, T], V) of the (complex) linear span of the collection of func¬ 
tions {xt -t G [0,T]}. By construction, relation (2.8) shows that we have an 
isometry between TCt and JCt, under which the correspondence Xt <—>■ xt 
holds true. Observe that under this isometry, we also have Mt <—> 1*, so, 
for every t G [0,r], the indicator function It belongs to JCt- It follows that 
JCt = LH[0,T],B[0,T],V). 

In the remainder of the paper we write and L‘^{[0,T],V) instead of 

L^(M,H(]R),//) and L‘^{[0,T],B[0,T],V), respectively. 


3. The transformation JCt — Ct- We have now associated three iso¬ 
metric Hilbert spaces with the fBm: the linear space TIt-, the frequency 
domain Ct and the space of integration kernels JCt- The aforementioned 
isometries between Ht and Ct and between TCt and JCt , determined by the 
relations Xt <—>■ et and Xt <—> xt-, respectively, induce a direct isometry 
between the function spaces JCt and Ct- We denote the map from JCt to Ct 
by U- 

Our first result gives an explicit analytic description of the isometry 
U :JCt —> Ct- The theorem states that it is a Fourier-type integral transfor¬ 
mation. The integration kernel is defined in terms of the function : M —> C, 
given by 

(3.1) ^(^) = /r(l-^)(j) ^7^0, 

U, ^ = 0. 

Here, as before, T is Euler’s gamma function and is the Bessel function 
of the first kind of order ia. Observe that, in fact, (p is defined on the whole 
complex plane. Moreover, property (2.3) of the Bessel function implies that 
p is an entire function. Evoking the well-known property 

(3.2) ^z"J,(z) = z"J,_i(z) 

of the Bessel function, one can easily see that the function p evaluated at 
Xt and the earlier introduced Fourier transform rht{X) are related by the 
identity 


(3.3) 


p{Xt) 


dmtjX) 

dmt{0)' 
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We need the following simple estimates for the function ip. The notation 
a < 6 means that a < cb, where c is positive constant that is universal or at 
least hxed throughout the paper. 


Lemma 3.1. For every A G M, the function u ^ p{uX) belongs to JCt 
and its norm satisfies 


(3.4) 


Il7^(-A)||y< 


1A|A|^-V2, 
1 V 


if <1/2, 

H>1I2. 


Proof. The fact that p is analytic implies that it is bounded in a 
neighborhood of 0. Using also that ^/zJu{z) is bounded for \z\ —> oo, we 
see that for real z, \p{z)\ is of order \z\^~^/‘^ for large \z\. So for H < 1/2 
we have \p{z)\ < 1 A \z\^~^/‘^, whence 

r\p{uX)\^dVu < A |A|2^-1) du < 1 A \X\^^-\ 

Jo Jo 

For H >1/2, it holds that \piz)\ < 1 V \z\^~^/‘^, so 

r\p{uX)\^dVu < V |A|2^-1) du < 1 V \X\^^-\ 

Jo Jo 

This completes the proof. □ 


The following theorem gives a complete description of the isometry U: JCt 
Ct- 


Theorem 3.2. The linear transformation U: JCt —> Ft is a Hilbert space 
isometry. For f G JCt, it holds that 

(3.5) Uf{X)= r fiu)p{uX)dV^ 

Jo 

for pL-almost all A G M, where p is given by (3.1). The class of functions 

(3.6) G Ct- J^\\p{-X)\\v\'iljiX)\p{dX) < ooj 
is dense in Ct and for G Cf- we have 

(3.7) U*fi{u) = U~^fi{u) = [ fi{X)p{uX)fi{dX) 

Jr 

for V-almost every u G [0,r]. Here U* denotes the adjoint of U. 
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Proof. By the Cauchy-Schwarz inequality and Lemma 3.1 


(3.8) 



for every A G M. Hence, the right-hand side of (3.5) defines a linear transfor¬ 
mation on ICt- Let us denote this transformation by A. 

Under the isometry ICt —> Ht, the indicator function It G ICt is mapped 
to the random variable Mt G Ht-, given by (2.5). Under the spectral isom¬ 
etry, Mt is mapped to the function A e-> rht{X) in Ct-, given by (2.2). So to 
prove (3.5), we have to verify that the mapping H, defined on ICt by the 
right-hand side of (3.5), coincides with the isometry U which is determined 
by the fact that Ult = rht for t G [0, T]. 

By (3.3), we have 



(3.9) 


So, indeed, Alt = Cht for every t G [0, T] and by linearity, U coincides with A 
on the set of simple functions in ICt- Now take an arbitrary / G ICt- The 
simple functions are dense in ICt, whence we can choose a sequence fn of 
simple functions such that /„ —> / in ICt- Then since U is an isometry, we 
have Afn = [//„—> U/ in Ct CL T^(^). On the other hand, (3.8) implies that 
Afn Af pointwise on M. But then Af and Uf must coincide for ^-almost 
all A G M, which proves (3.5). 

Both the isometry ICt —> Ht and the spectral isometry Ht —> Ct preserve 
inner products, so the same holds for their composition U '.ICt —> Ct- This 
implies that U is unitary, that is, that U~^ = U*, where U* : Ct —> ICt is the 
adjoint of U, determined by the relation {Uf, = (/, U*'ip)v for all / G ICt 
and V’ £ Tt- Using (3.5), we see that, for ip G Ct, we have 


{Uf,iP)^= [ UfiX)iP{X)fi{dX) 
Jr 



For Ip G C'rp, we may interchange the integrals, since by the preceding lemma 
and Cauchy-Schwarz, 



which is finite by definition of C'j-. It follows that 



which proves (3.7). 







SPECTRAL THEORY FOR THE FBM 


9 


It remains to prove that is dense in Ct- Let the S be the Schwarz 
space of rapidly decreasing functions on M, that is, C^-functions / on M such 
that for all m,n, the derivative satisfies \x |m|j(n)^^)| ^ Q ag |3;| ^ 
and let St be the space of Schwarz functions with support in [0,r]. By the 
preceding lemma, it clearly holds that S H Ct C so it is enough to show 
that S n Ct is dense in Ct ■ 

Fix t e [0,r] and choose a sequence Wn of C°° probability densities such 
that supp(tc„) C [0,t] and such that the associated probability measures 
converges weakly to the dirac measure 6t concentrated at t. Dehne 





Wnis)ds, 


0 , 


u<t, 
u G {t,T]. 


Then /„ is a C'^-function with compact support, so /„ € St- For fixed 
A G M, the function s i—> l<i(A) is bounded and continuous, whence the weak 
convergence implies that 


fnW= [ is{X)Wn{s)ds^ [ is{X)St{ds) = it{X) 
Jo Jo 

for every A G M. Observe also that, for A G M, 


l/n(A)|^< 



is(A)|^t(;n(s)ds < 1 A 


1 

A 2 - 


By dominated convergence, it follows that /„ —> It in Since the 

Fourier transform maps St into S, the functions fn belong to 5 H Ct- So for 
every t G [0,T], It is the L^(/x)-limit of a sequence of functions in 5 n Ct- 
Since Ct is the closure in of the linear span of the functions If, 

t G [0,T], this shows that S H Ct is, indeed, dense in Ct- □ 


Observe that since we have 


(3.10) Jy^iz) 



•^- 1/2 (^) 



Zy^O 


and r(l/2) = it holds that (p{z) = in the standard Brownian motion 
case H = 1/2. So in that case, the map U is simply the Fourier transform. 
For general H G (0,1), we can view it as a fractional version of the Fourier 
transform. 

It seems worth mentioning that fractional integration theory enters in the 
present context via the simple observation that the Fourier transform rht 
of the kernel ruf, defined by (2.1), is expressible in terms of the fractional 
integral of order 3/2 — H oi the function exp(iuA). This can be 
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seen by comparing (2.2) with formula 9.1.10 of Samko, Kilbas and Marichev 
(1993). Specifically, we have that 


mt{X) = 


1 




r{H + i/2) °+ 

By (3.3), it follows that, for the Fourier-kernel of the map U, it holds that 




1 




r(i? + i/2) °+ 

Hence, using fractional integration by parts, we see that, for / G ICt, 

1 ‘■T 


Uf{X) = 


T{H+ 1/2) 

1 

T{H+ 1/2) 


Jo 

.F(u'/2-^4/_2-^/(u))(A), 


provided, of course, that the fractional integral of order 1/2 — H of / and the 
Fourier transform (denoted by T) exist. The composition rule of fractional 
integration operators implies that, for -0 G Ct, 

U-^'i/it) = T{H + l/2)4_r^/^(u^-i/2^-V('«))(i)- 


Note, for instance, that for a deterministic integrand / €It, the latter ex¬ 
pression for U~^, in combination with the spectral isometry, yields the re¬ 
lation 

r f{u) dXu = T{H + 1/2) r I^:^^\u^-^/^f{u)){t) dMt, 

Jo Jo 

where M is the fundamental martingale. For / = 1^, this reduces to the mov¬ 
ing average representation (2.7). In general, the expression of the operators 
U and in terms of Riemann-Liouville operators, can be very useful 

for the evaluation of the transforms in concrete cases, since many explicit 
formulas for fractional integrals are known. We will, however, not need this 
connection in the present paper. The proofs of our results do not use any 
fractional calculus. 

Relation (3.5) gives an analytic description of the functions in Ct- In 
particular, it allows us to prove that every function in Ct is the restriction 
to M of an entire function. Strictly speaking, the elements of Ct are, of 
course, equivalence classes of functions. Two functions represent the same 
element if they coincide /r-almost everywhere. Theorem 3.2 implies that 
every equivalence class can be represented by an entire function. 


Corollary 3.3. Every element of Ct has a version that is the restric¬ 
tion to ffi of an entire function. 
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Proof. For / £ JCt, consider the complex function 

rT 


(3.11) 


/o 


f{u)<f{uz)dVu. 


Since (p is entire, this function is well defined and easily seen to be continuous 
on C. To prove that the function is analytic, consider a closed path 7 in the 
complex plane. By Fubini’s theorem and Cauchy’s theorem, 

. ' ^ r ^ [ . , 

/ \Io Jo (/ ^ 

Hence, by Morera’s theorem, the function defined by (3.11) is entire. □ 

In the remainder of the paper, if we consider an element ip £ Ct, we will 
always assume this to be the smooth version. 


4. The reproducing kernel on Ct- As was established in Theorem 3.2, 
the transformation U is of Fourier-type, generated by the Fourier kernel p. 
This motivates us to introduce the function St on M x M, defined by 


fT _ 

5r(a;,A)= / ip{uLo)ip{uX) dVu, 
Jo 


where (p is, of course, given by (3.1) again, and V by (2.6). The Cauchy- 
Schwarz inequality and Lemma 3.1 imply that St is well defined. Moreover, 
(3.5) implies that, for fixed a; £ R, the function A 1 —> St^, A) is the image un¬ 
der U of the function u 1 —> p{uuj). In particular, we see that A 1 —> St^, A) be¬ 
longs to Ct for every a; £ R. It clearly holds that 5 t ( w , A) = St{X,u>). Since 
99 ( 0 ) = 1 , relation (3.9) implies that S'r(0,0) = Vr and S't(0, A) = mT(A). 

The following theorem states that St acts as a reproducing kernel on the 
spectral space Ct, turning it into an RKHS. 


Theorem 4.1. For every ■0 £ Ct, we have 


iIj{\)St{uj, X)p{dX) = V'(a;), 


for all uj ^ . 


Proof. Suppose first that ?/) £ 5 n Ct, where S is the Schwarz space of 
rapidly decreasing functions. Then by Fubini’s theorem and Theorem 3.2, 

[ V'(A)S't(u;, A)/r((iA) = [ iIj{X)( [ p{uijj)p{uX) dvi\p{dX) 

Jm. Jk. \Jo ) 

= p{uuj) '4!{X)p{uX)p{dX)^ dVu 
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= [ ip{uuj)U ^'il){u)dVu 
Jo 


= [/([/-= 


for ;U-almost all u; G M. The interchanging of the integration order is justihed 
by the fact that -0 is rapidly decreasing. 

Now let G Ct be arbitrary. Since S n Ct is dense in Ct (see the proof 
of Theorem 3.2), we can choose functions G 5 H Ct such that 'tpn ^ iIj 
in L^(//). By the remarks preceding the theorem, the function A e-> St{u:, A) 
belongs to for fixed a; G M. So by Cauchy-Schwarz, we have 


V^n(A)5T(u;,A)//(dA)- / ^|JiX)ST{co,X)^i{dX) 


< 


— Yn 


,||5r(w, OIU ^ 0 


for every cu G M. By the preceding paragraph, the first integral on the left- 
hand side equals for ^-almost every oj. So the functions V’n converge 

^-almost everywhere to the function 


UJ 


V’(A)5T(cu,A)//(dA), 


and they converge in L?‘{^Cj to ij:. Then the two limits must coincide /r-almost 
everywhere. Since both functions are continuous (Corollary 3.3), the proof 
is complete. □ 


A first simple consequence of the RKHS structure is that the “kernel 
functions” X^ St{u}, X) span Ct- 

Corollary 4.2. The space Ct is the closure in of the linear span 

of the collection of functions {A e-> ^^(w, A): a; G M}. 

Proof. We already noted in the beginning of the section that every 
function ^^(w,-) belongs to Ct, so the closure of the linear span of{Ai—> 
S'T(ti;, A): u; G M} is certainly contained in Ct- 

To prove the inclusion in the other direction, take G Ct and suppose 
that if} is orthogonal to every kernel function 5r(w,-), that is, 

/ V'(A)5r(w,A)/x((iA) = 0 

Jk. 

for all a; G M. Then, by the reproducing property of St, we see that if 
vanishes //-almost everywhere. □ 
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It seems useful to briefly discuss the relation between the preceding frequency- 
domain results and the so-called “time-domain RKHS.” The latter space is 
constructed by associating to every element H G Tij- a function 1 1 —> KHXt 
on [0,r]. These functions are the elements of the time-domain RKHS and 
the inner product of two functions 1 1 —> EHXt and 1 1 —> EH'Xt is defined 
as E,HH'. By construction, the resulting Hilbert space is isometric to TCt, 
and the covariance function r{s, t) = is the reproducing kernel on the 

space. 

The following theorem clarifies the relation between the two reproducing 
kernel Hilbert space structures. 

Theorem 4.3. The time-domain RKHS is the closure of the linear span 
of the collection of functions {t lt(A): A G M} with respect to the inner 
product 

{t 1-^ 11 —> it(A)) = A). 

Proof. The spectral isometry shows that the time-domain RKHS is 
given by all functions 1 1 —> {fj, it)^ on [0,r], where if runs through Ct, the 
inner product of two elements t ^ and 1 1 —> (^,1^)^ being given by 

(V’jOai- Corollary 4.2, it follows that the time-domain RKHS is the clo¬ 
sure of the linear span of collection of functions {t i—> (S't(A, : A G M}. 

By the reproducing property, we have (^^(A, •), = lt(A) and the inner 

product between the functions 1 it(A) and 1 1 —> it(w) equals St{oj, A). □ 

We note that the moving average representation of the fBm implies that 
the inner products in the time-domain RKHS can be expressed in terms of 
Riemann-Liouville fractional integration operators [cf., e.g., Hult (2003)]. 
Theorem 4.3 thus yields an expression for the reproducing kernel St of Ct 
in terms of fractional integrals. In the present paper this connection plays 
no further role, however, and we will now return to the frequency domain 
analysis. 

5. Differential equations. Our next intension is to obtain an explicit 
analytic expression for the reproducing kernel St- This goal is achieved 
in Corollary 6.2 below. In this section we will show that the Fourier ker¬ 
nel (p is subject to a “recursion relationship” that is also encountered in 
the study of other types of processes with stationary increments [see, e.g., 
Kailath, Vieira and Morf (1978), where a short account can be found of 
the classical result of Krein (1955) concerning the “signal plus white noise” 
model]. 
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The Bessel function satisfies the second order ordinary differential 
equation 

(5.1) Jui^) + ~ = 0- 

In view of the representation (2.2), this yields the following differential equa¬ 
tion for the functions rht € Ct- 


Lemma 5.1. For every A £ M, we have 


d‘^mt{X) _ dmt{X) 
dt‘^ dt 




dmt{X) 

dt 



where 'yt = {H — lj2)jt. 


Proof. By (2.2), we have rht{X) = cxfu{z)Ju{z), where cx is a constant 
(depending on A), z = At/2, v = 1 — H and fuiz) = z'^ ex.p(iz). It is easily 
verified that 

d t t \ _ fu{z){v + iz) ^ ^ _ fu{z){{v + izf - u) 

’ t ’ ’ t2 

The differential equation is now a straightforward consequence of (5.1). □ 


This differential equation for the functions 1rht{X) gives rise to a dif¬ 
ferential equation for the Fourier kernel <f{tX). We present this as a system 
of equations for the function P{t, A), defined for t > 0 by 

(5.2) P{t,X)=ip{tX)^JdVt/dt, 
and its reciprocal 

(5.3) P*{t,X) = e^^*P{t,X). 

Note that 

(5.4) 5 ^( 0 ;, A) = [ P{t, X)P{t, X) dt. 

Jo 

We also observe that since the functions rht correspond to the fundamental 
martingale M under the spectral isometry, the functions 


f P{s,-)ds 
Jo 

correspond to an ordinary Brownian motion. In particular, we have that 

[ f P{u,X)du [ P{v, X) dv fi{dX) = s At. 

Jr Jo Jo 
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If we interchange the integrals on the left-hand side and differentiate, we get 
the formal expression 

/ P{s,X)P{t,\)fj,{dX) = S{t — s). 

Jr 


In this sense, the functions P{t, ■) are orthogonal with respect to the spectral 
measure 

The following theorem shows that the orthogonal functions satisfy Krein’s 
continuous analogue of the usual recurrence formulas for orthogonal poly¬ 
nomials on the unit circle. 


Theorem 5.2. For every A G R, the function P{t,X) and its reciprocal 
P*{t,X), defined by (5.2) and (5.3), satisfy the equations 

dP{t, X) 


dt 


= iXPit,X)-jtP*it,X) 


and 


dP*{t,X) 

Ft 


= -ItPitA), 


where ^t = {H — l/2)/t. 


Proof. First of all, let us express P{t,-) and P*{t,-) in terms of the 
function rht given by (2.2). We have drhfiX) = ip{tX) dVt and Vt = 

[see (3.3) and (2.6)], hence, 

t^~^Pm'fiX) 


(5.5) 


P{t,X) = 


^{2-2H)dl 

where the prime denotes differentiation with respect to t. It follows that 

t^-V 2 (m')*(A) 


(5.6) 


P*it,X) = 


^{2-2H)d]j ’ 


where (m()*(A) = exp{iXt)m[{X) is the reciprocal of m[. Observe that 

- y^i(A) = e*"*/2|(e-*"*/2mt(A)), 


^m((A) — (A) = e*^*e*^*/2^(e“*^*/2mf(A)) 

= e*^'/2^(e-*^*/2^t(A)) 

= mt(A) - ymt(A). 


so that 











16 


K. DZHAPARIDZE AND H. VAN ZANTEN 


Since rht is self-reciprocal, that is, ml{X) = exp(iXt)mt{X) = mt{X), the 
latter identity implies that {m[)*{X) = rh^{X) — iXrht{X). Combining this 
with (5.6), we find that 

(5.7) o.>. 

^(2-2H)cP„ 

The first statement of the theorem now follows from differentiation of (5.5), 
taking Lemma 5.1 and (5.5) and (5.7) into account. Similarly, the second 
statement is obtained by differentiating (5.7). □ 

6. ChristofFel-Darboux formula. As in the theory of orthogonal polyno¬ 
mials and their continuous analogous, the “recurrence relations” presented 
in Theorem 5.2 allow us to derive a closed-form expression for the repro¬ 
ducing kernel St [cf., e.g., Grenander and Szego (1958), Section 2.3, and 
Kailath, Vieira and Morf (1978), formula (48)]. 

Theorem 6.1. Let P{t,X) and its reciprocal P*{t,X) be defined by (5.2) 
and (5.3) . For all T > 0 and w, A e M, we have 

(6.1) i{X-uj)ST{Lo,X) = P{T,lo)P{T,X) - P*{T,uj)P*{T,X). 

Proof. We view the left-hand side and right-hand side of (6.1) as func¬ 
tions in T. Recall that we have (5.4). Using Theorem 5.2, a straightforward 
calculation shows that 

i{X - eo)P{^)P{fi X) = ^( P(^P(t, A) - PH^P*{L A)). 

This shows that the functions on both sides of (6.1) have the same derivative 
with respect to T, which implies that their difference is independent of T. 
So for every T > 0, we have 

(6.2) i{X - oj)St{uj, A) = P{T,u})P{T, A) - P*{T,oj)P*{T, X) + C(w, A) 

for some constant C{uj,X), and it remains to show that C{uj,X) = 0. 

For H < 1/2, the functions t e-> P{t, A) are bounded for every A G M and it 
holds that P{0, A) = 0. So in this case we can let T —> 0 in (6.2) to see that 
C{lo, A) = 0. Since the Bessel function J^, is analytic in n [see Watson (1944), 
page 44], (6.2) shows that as a function in H, C{io, A) can be extended to an 
analytic function on the open disc of diameter 1, centered at 1/2. We just 
saw that it vanishes for all H in the interval (0,1/2), whence a standard 
result from complex function theory implies that it vanishes on the entire 
disc [cf. Rudin (1987), Theorem 10.19]. In particular, C{lo,X) = 0 for all 
FT G (0,1). □ 
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In combination with the explicit expression that we have for the orthog¬ 
onal functions the preceding theorem yields an explicit analytic ex¬ 

pression for the reproducing kernel St- 


Corollary 6.2. 
resentation: 

(i) For uj 

St{2lo,2X) 

5t(0,0) 


The reproducing kernel St admits the following rep- 


J.h{Tw)Ji-h{TX) - Ji.h{Tw)J-h{TX) 
T{X-io) 


(ii) For u) G 
St{2lo,2uj) _ 
5t(0,0) 


/ (yTT _ -I \ 

X [jI_h{Tuj) + -j^J-h{Tw)Ji.h{Tw) + JIh{Tuj)] • 


Proof. Part (i) follows by straightforward calculations from the pre¬ 
ceding theorem, the definition (5.2) of P(t,X), the explicit expression (3.1) 
for If and the fact that dVt/dt = (2 — 2H)Vt/t. 

To prove part (ii), we note that St is analytic and, in particular, contin¬ 
uous, so we may derive an expression for St{co,Lo) by letting A —> a; in the 
expression that we found in part (i). It suffices to observe that as A —> w, we 
have 

J-H{tio)Ji-H{tX) — Ji-H{tco)J-H{tX) 

X — CO 


= J-H[tCO) --- Ji-H[tCO) - 


X — CO 

d d 

J-H{tCo) — Jl-H{tco) — Ji-H{tco) — J-H{tuo) 


X — CO 


/ \ 

= t f Ji_jj{tLo) H-- J—Hitu!)Ji—H{tUj) j‘^jj{tco) 


too 


In the last step we have used the recurrence formulae 

(6.3) 
and 

2v + 2 

(6.4) Jy+2{z) = - Ju+i{z) - Ju{z) 


■^Ju{z) = -Ju{z) - Jv+i{z) 
dz z 


[see Watson (1944), page 45]. □ 


2 ; 
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7. Orthonormal basis in Ct- By the reproducing property, the inner 
product of the kernel functions A i—> 5 t(2w, A) and A i—> 5 ^( 20 ;', A) in Ct is 
given by S't(2u;, 2a;^). Hence, by Corollary 6.2, these functions are orthogonal 
in Ct if Tw and Tu>' are different zeros of Ji-h- In this section we prove 
that if we let Tco range over all zeros of Ji-h, we obtain an orthogonal basis 
of Ct- 

We first recall some facts about the zeros of the Bessel function of the first 
kind [see, e.g., Erdelyi, Magnus, Oberhettinger and Tricomi (1953), Section 7.9]. 
For 1 / > - 1 , the Bessel function Ji, has a countable number of positive 
zeros that can be ordered according to magnitude. We denote them by 
Ai /,1 < Xu ,2 < • • • . For positive u, the function Ju satisfies Ju{0) = 0 and its 
negative zeros are given by —Xu,i > —Xu ,2 > ■ ■ ■ ■ Hence, for > 0, the ze¬ 
ros of Ju can be ordered as • • • < Xu,-i < Xu,o = 0 < Xu,i < ■ ■ To prove the 

completeness of the system of orthogonal functions, we need the following 
lemma, which is a consequence of the multiplicative decomposition of the 
Bessel function, or, more precisely, of the formula 

T (7^ ~ ^ >2 _ ^2’ 

see Erdelyi, Magnus, Oberhettinger and Tricomi [(1953), formula 7.9.3]. 


Lemma 7.1. Let Ju be a Bessel function of the first kind of nonnegative 
order iy>0 and let ■ ■ ■ < A_i < Aq = 0 < Ai < • • • be its real zeros. Then 


Ju+l (- 2 ) 
Ju{z) 


Ju+l{w) 

Ju{w) 



Z — W 
■2') (An 


w) 


for all 2 ;,u; G M, for which the expressions are well defined. 


Proof. Since the zeros of Ju satisfy A_n = —An for all n G N, it holds 
that 


E 

n^O 


Z — W 


(An z){Xn w) 


Z — W 


+ 


z — w 


(An . 2 ) (An w) (An T 2;) (An T w') 


Now observe that we have the identity 

a — b a — b 

+ ■ 


2a 


2b 


{c — a){c — b) {c + a){c + b) — a? c^ — b'^' 
for all a, 6, c G K for which the expression makes sense. It follows that 
z — w 2z 2w 


E 

n^O 


= E 


{Xn-z){Xn-w) ^A2-2:2 


-E 
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By (7.1), the right-hand side is equal to 

Ju+i{z) _ Ju+i{w) 

Ju{z) Jv{w) 

and the proof of the lemma is complete. □ 


We can now present the orthonormal basis in Ct and the associated ex¬ 
pansion formula. 


Theorem 7.2. Let ■ < u-i < wq = 0 < wi < • • • be the real zeros of Ji-h 
and, for n G Z, define the function ipn on M by 


and put 


i’nW 


ST{2u;n/T,X) 

||5T(2n;„/r,.)IU’ 


(7.2) 


[Vt, 


7^ 0 ? 

iOfi — 0. 


The functions ipn form an orthonormal basis of Ct ond every function xjj E 
Ct can be expanded as 

nez V / 


the convergence taking place in L^(/r). 


Proof. By the remarks in the beginning of the section, the functions ifn 
are orthogonal, and they have unit norm by construction. Let us prove that 
the system is complete. 

By Corollary 4.2, it suffices to show that every kernel function A 
ST{aj,X) is in the closure of the linear span of the tpn- We claim that, for 
cj, A G M, 

S'r(2a;,2A) = ^ ipni2uj)fJni2X). 

n^'L 


To prove this, note that we have 


(7.3) 


^ V’n(2w)V’n(2A) 


5T(2cj,0)5r(0,2A) 

5t(0,0) 


5'j'(2u;, 2w„/T)5'7'(2u;^/T, 2A) 

ST{2u:nlT,2unlT) 
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By Corollary 6.2, the sum on the right-hand side equals 


^^ Ji-h{T\)Ji_h{Tlo) 


JT{X-uj) 


E 

nj^O 


{LOn-TX){u;n-T(jj)' 


where c = (2 — 2H)T^{1 — H)Vt- By the preceding lemma, the sum in the 
last display equals 

1 / J2-HiT\) J2 -h{Tuj)\ 

T{\-oj)\Ji.h{T\) ~ Ji-h{Tuj))' 

In view of (6.4), multiplication by Ji-h{T\)Ji-h{Tuj) yields 

J2-h{T\)Ji.h{Tu:) - J2-h{Tu:)Ji-h{T\) 

T{X-uj) 

_ J.h{Tu2)Ji.h{TX) - J.h{TX)Ji.h{Tuj) 

T{X-uj) 

9 _2/7^ 

Hence, using Corollary 6.2 again, we see that the sum on the right-hand side 
of (7.3) equals 

St{‘2uJ,2X) — ^ gtr(A-a;) Ji_h{TX)Ji-h{Tlo). 

Therefore, our claim follows from the fact that the second term in the last 
display is equal to 

S’'7’(2u;, 0)S''7’(0,2A) 'rhj'{2X)fhT{‘2uj) 

S'r(0,0) mr(0) 

[recall that £' 7 ’( 0 ,A) = mT{X) and 5 t(0,0) = Vt = m7’(0)]. To check this, 
evoke expression (2.2). 

So, indeed, the functions form a complete, orthonormal system. It 
follows that every '0 G Ct can be written as ijj = '0n)^V'n- By the re¬ 

producing property. 




J^^l;iX)ST{2uJn/T,X)^i{dX) _ i;{2oJnlT) 


||5T(2u;„/r,.)IU ||5T(2a;„/r,.)IU' 

Another application of the reproducing property yields 

The explicit expression in (7.2) of the normalizing factor a{ujn) follows from 
the second assertion of Corollary 6.2. □ 

























SPECTRAL THEORY FOR THE FBM 


21 


We remark that, instead of the zeros of Ji-h, we can also use the zeros 
of J-H to obtain a second orthonormal basis of Ct- Since the reasoning is 
completely analogous to the Ji-h case, we omit the details and mention only 
that all the consequent results of this paper can be easily reformulated in 
terms of the zeros of J-h [like in the case H = 1/2 oi the Brownian motion, 
where there exist expansions in terms of the zeros of the sine and the cosine, 
cf., e.g., Yaglom (1987), Section 26.1], 

Using the isometry U: JCt —> Cti we can now easily obtain an orthonormal 
basis of the function space ICt = L‘^{[0,T], V) and the corresponding series 
expansion. Note that in the case H = 1/2 it reduces to the Fourier series 
expansion, given for / G L^[0,r] by 


Corollary 7.3. Let • • • < oj-i < cuq = 0 < cui < • • • he the real zeros 
of Ji-H o,nd let (p he given hy (3.1). Then the functions 


u t 


(p{{2uJn/T)u) 


n G Z, 


||^((2o../r).)||u’ 

form an orthonormal basis of ICt = Lp‘{\fl,T],V). Every f & ICt can he ex¬ 
panded as 


f{u) = Y. 


u f{2u)n/T)ip{{/2uin/T)u) 
M{2u:n/T)-)\\l 


the convergence taking place in L^([0,T],U). 

Proof. Just note t hat St {(jJ, •) is the image under the isometry U: ICt — > 
Ct of the function u i—> p{uu>). □ 

8. Paley-Wiener expansion. In this section we use the orthonormal basis 
of Ct to obtain a series expansion of the fBm. Paley and Wiener (1934) use 
a series expansion to introduce the standard Brownian motion, which they 
call the “fundamental random function” [see also Levy (1965), Section 13 of 
the Complement]. They first consider the series 

tG[ 0 , 27 r], 

where the are i.i.d., complex-valued, standard Gaussian random vari¬ 
ables. This series corresponds to white noise, but the series does not converge 
in the usual sense. So instead they consider its formal integral 

^int _ 


T- 

nEZ 


in 


-Zn- 
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The latter series is shown to converge almost surely and is taken as the 
definition for the (complex-valued) Brownian motion. 

We can now present the extension to the case H ^1/2 of this expansion. 
It is natural to consider a complex-valued fBm. This is a centered, complex¬ 
valued Gaussian process X with covariance structure 

EXsXt = - |s - 


Theorem 8.1. Let ■ ■ ■ < a;_i < wq = 0 < wi < • • • be the real zeros of Ji-h 
and let Z^, n G Z, he independent, complex-valued Gaussian random vari¬ 
ables with mean zero and variance E|Z„|^ =cj^(a;„), where a‘^{iOn) is given 
by (7.2). Then, with probability one, the series 

2iu}ntlT _ 1 

E 9- ,T 

converges uniformly in t € [0,T] and defines a complex-valued fBm with 
Hurst index H. 

Proof. By Theorem 7.2, we have 

n& \ J- / 

where ijjn is a complete orthonormal system in Ct- It follows that, for G 

[o,U, 

(8.1) (i., i,)„ = (^) i,(^), 

By the spectral representation of the fBm, the left-hand side of the display 
equals the covariance function of the fBm. Hence, the equality shows that 
the series in the statement of the theorem converges in mean square sense 
for every t G [0,r] and defines a Gaussian process with the same covariance 
structure as the fBm. 

The fact that the series converges uniformly with probability one can be 
deduced from the Ito-Nisio theorem [cf., e.g., Ledoux and Talagrand (1991)]. 
See, for instance, the proof of Theorem 4.5 of Dzhaparidze and Van Zanten 
(2004) for details. □ 

For H = 1/2, we have ujn = mr [see (3.10)] and = 1/T, so, indeed, 

we recover the classical Paley-Wiener expansion of the ordinary Brownian 
motion in this case. 
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The Paley-Wiener theorem for the fBm shows that the (complex-valued) 
fBm on [0,1] can be viewed as the formal integral of the process 

In view of (7.2) and the fact that z y/zJu{z) is bounded, the latter can 
be seen as a random signal in which the weight of the component with 
frequency ojn is (up to a constant) approximately equal to ^. 

The real-valued version of the expansion is as follows. 


Corollary 8.2. Let X, {Yn)ne^ o,nd {Zn)n&i be independent, real¬ 
valued Gaussian random variables with mean zero and variance 

= —, EK^ = 

Vt 2 

where (T^(a;„) is given by (7.2). Then, with probability one, the series 


tx + E 

n=l 


sin 2u}nt/T 
^^nlT 


u + E 

n=l 


{cos 2uJnt/T — 1) 
^nlT 


Zr, 


converges uniformly in t € [0,T] and defines a real-valued fBm with Hurst 
index H. 


Proof. Note that the terms indexed by n and —n in the sum in equa¬ 
tion (8.1) are complex conjugates, hence, 

{i„ i,)„ = + 2 f; ^"K)s(i, (^) i, (^) ). 

Since 



s\\l2u}nS/T SYn2iu)nt/T {cOs2uJnS/T — 1) {cOs2iOnt/T — 1) 
2WnlT 2lWnlT 2.lWnlT 2uJnlT 

this shows that the series in the statement of the corollary converges in mean 
square sense and that the resulting process has the same covariance structure 
as the fBm. Uniform convergence with probability one can be argued as in 
the proof of the theorem. □ 


In the paper Dzhaparidze and Van Zanten (2005) we proved that a related 
series expansion of the fBm is rate-optimal in the sense of Kiihn and Linde 
(2002). This means that the rate of uniform convergence is N~^^J\ogN. 
Exactly the same reasoning as in Dzhaparidze and Van Zanten (2005) can be 
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used to prove the optimality of the Paley-Wiener expansion of Theorem 8.1. 
The main idea is simply to use the asymptotic properties ~ nvr and 
for n ^ oo [cf. Watson (1944)] to estimate the size of the 
terms in the expansion. We omit the details and just give the optimality 
result. 


Theorem 8.3. The expansion ofthefBm of Theorem 8.1 is rate-optimal. 
It holds that 


E sup 

ie[0,T] 


E 

jn|>Af 


‘liixinlT 





Let us mention that related optimal series expansions for the so-called odd 
and even parts of the fBm, and, consequently, also of the fBm itself, can be 
found in Dzhaparidze and Van Zanten (2004). Compared with the Paley- 
Wiener expansion of Theorem 8.1, the representation of Dzhaparidze and Van Zanten 
(2004) has the drawback that it requires both the positive zeros of J\-h and 
of J-u- 

Another expansion optimal in the sense of Kiihn and Linde (2002) is ob¬ 
tained by Ayache and Taqqu (2005). Their construction is of a completely 
different type. It involves a doubly indexed array of i.i.d. Gaussian random 
variables weighted by functions 1 1 —> nifldt — k) — '!'//(—A:)), de¬ 

fined in terms of the Fourier transform of an appropriate mother wavelet T. 

Finally, we want to mention the possibility to extend the expansion re¬ 
sults to the fractional Brownian sheet (fBs). This can be achieved by taking 
suitable tensor products like in the paper of Dzhaparidze and Van Zanten 
(2005), where our earlier double series expansion [cf. Dzhaparidze and Van Zanten 
(2004)1 is extended to the fBs. On extending the Paley-Wiener expansion of 
Theorem 8.1, the construction is analogous and even simpler because we now 
have only one sequence of Bessel zeros. The resulting expansion of the fBs 
is again rate-optimal, as can be shown by using the asymptotic properties 
of the Bessel function and its positive zeros. 


Acknowledgment. Thanks to Michael Lifshits for pointing out some silly 
miscalculations. 
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